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Abstract 

We present functional Schrodinger representations of Holstein-Primakoff bo- 
son and slave-boson theories for the Heisenberg Hamiltonian and the t — J 
Hamiltonian respectively. Based on these representations we obtain the dis- 
persion relations of magnons for two dimensional antiferromagnets. By ap- 
plying the functional Schrodinger representation of the Holstein-Primakoff 
boson theory to the Heisenberg Hamiltonian, the exchange energy is correctly 
predicted and the self-energy of quasi-hole is obtained. From the use of the 
functional Schrodinger representation of the t — J Hamiltonian it is shown 
that at half-filling the dispersion relation obtained from the slave-boson the- 
ory leads to that obtained from the Holstein-Primakoff boson approach. 

PACS numbers: 71.10.Fd, 71.27.+a, 75.30.Ds 
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I. INTRODUCTION 



Although widely used in the area of high energy physics, the functional Schrodinger 



number of application to condensed matter physics appeared in the literature .113113 Earlier 
we treated the t — J Hamiltonian by introducing a slave-boson approach to the functional 
Schrodinger picture(FSP) theoryB for the two-dimensional systems of antiferromagnetically 
correlated systemsEl. The Holstein-Primakoff boson representation is often used to describe 
the broken symmetry phases of the quantum Heisenberg antiferromagnet. Lately Changjl! 
introduced a generalized Holstein-Primakoff representation of the t — J model Hamiltonian 
in order to describe a higher order(second order) effect on spin waves and antiferromagnetic 
spin polarons by allowing a systematic perturbative expansion. Here we present the func- 
tional Schrodinger representations of both the Holstein-Primakoff@ boson and slave-boson0 
theories for the Heisenberg Hamiltonian and the t — J Hamiltonian respectively. Based on 
these representations we derive magnon dispersion relations and present comparison between 
the two approaches. 

II. HOLSTEIN-PRIMAKOFF BOSON THEORY OF HEISENBERG 
HAMILTONIAN BY FUNCTIONAL SCHRODINGER REPRESENTATION; 

ANTIFERROMAGNETIC MAGNON 

To allow symbol definitions for later use, here we choose a brief review on a generalized 
approach of perturbatively treated Heisenberg Hamiltonian based on a Holstein-Primakoff 
transformation!^. The Hilbert space of present interest is spanned by the states, |oi) ® \a 2 ) ® 
•• •<£) \cij) Cg)- • -® \ajsf), where |oj) G {c^ |0)j), cjj0)j, |0)j} with Cj a \0)j = 0. The vacuum state 
is given by |0) = <8>|0)j which satisfies c Jcr |0) = 0. j is the site index and a, the spin index. 
N is the total number of lattice sites. With the local Hubbard operators Xj b = \a,j){bj\, the 
Heisenberg Hamiltonian for the two-dimensional system of antiferromagnetically correlated 
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electrons is written) 



H, 



J 



E( \ r <J—(T \r — a(7 \r(j(j \r — a — <j 

^Aj A 3 - -Aj A j 



(1) 



where J is the Heisenberg coupling constant, (ij) stands for summation only over nearest 
neighbours. The above Hamiltonian includes the contribution of the jfiiUj term that ap- 
pears in the usual t — J Hamiltonian, where n, = J2a c la c i(r- Introducing two commuting 
boson(magnon) operators (cij,bj) and an anticommuting fermion(hole) operator fj for a bi- 
partite lattice made of sublattices, A and B, two sets of Holstein-Primakoff representations^ 
for the local Hubbard operators are given in the table below. 



Hubbard operator 



oo 



X 



xf 



X 



to 



X 



oi 



X 



X 



TT 



X 



X 



A sublattice 



fjfi 



H \ • - i («>; • /;/,•) 



i 



2s wl — 7T- I aia,- 



2s 



2s V J ■? 



f]fj)f: 



2s \ a ] a j + fjfj 



2sa! 



1 _ _ I a T a . + /Tj. 



_B sublattice 



2s 



/2s 



For sublattice A, a,- (a!-) is the annihilation(creation) operator of spin-up boson (magnon 



excitation) and fj(fj), the annihilation (creation) operator of spinless fermion (hole ex- 
citation) at site j. s is the spin quantum number. Likewise, for sublattice B, bj(bj) is 
the annihilation(creation) operator of spin-down boson(magnon excitation) and fj(fj), the 
annihilation(creation) operator of spinless fermion(hole excitation) at site j. 

Using the local Hubbard operatorsEl in the table above, we derive from Eq. ([!]), 

Hj = j£ jsatyi - 1 (a\a t + fj^l- 1 (bfa + fj f 3 



-~ alaitfjbj + sJ 1 - — (a\ai + fjfi)aj 1 - - (fetfy + fjfj)bj 



-2s' 



(2) 



The transformed Heisenberg Hamiltonian above is perturbatively treated by using the fol- 
lowing expansions; 



Wi-iO^ + ZM) = 1 - h + f}fi) + {H-O-) 



(3) 



2s V •' 4s 

We consider terms up to order - in Eq. for the evaluation of Eq. (Q). By first taking 
the Fourier transformation of Eq. (0) and then the Bogoliubov transformation, we obtain 
the following four terms; 1) the free Hamiltonian Hq, 



.<4 



k 

-Js 2 zN. 



i-Iiv, 

2s N ^ \Jsz 



- 1 



c k c k + Ad*) + Jsz^ 



^k 



1 - 2/i/ k 



(4) 



with AT, the total number of spins. 2) The magnon-magnon interaction, H 



M-M, 



H 



M-M 



Jz 2 v 

ki,k 2 



^k^k 2 



(Jsz) 2 

(J sz \t ) C ici C ki ct k 2 Ct k 2 



( c ki c ki c k 2 c k 2 + ^ki ^ki G^k 2 <^k 2 



+2 (l + ^Wc kl dU 



(5) 



3) the magnon-hole interaction, Hm-h, 



H 



M-H 



Jz 2 



a; 

E j^ftj* ( c k 2 c k2 + 4 2 ^k 2 ) 



ki,ki 



(6) 



and 4)the hole- hole interaction, H H _ H , 



Jz 2 



Hh-H — -rz E Tq/ki-q/ki/^+q/l 



2 iV. , 

ki,k 2 ,q 

where cj£ = JszJl - 7^, 7 k = -jrE^' 5 



k 2 , 



(7) 



7_k with the number of nearest neighbors. 
For a square lattice we have 7^ = ~ [cos(k x a) + cos(k y a)]. The last two terms in Eq. are 
contributed from the 2 n i n j term that appears in the original t — J Hamiltonian. This term 
was not considered in the work of ChangEl. Ck(c k ) and rfk(^L) are the annihilation(creation) 



operators of the chargeless bosons(quasimagnons) of momentum k, in association with two 
sublattices A and B respectively. /k(/k) is the annihilation(creation) operator of a spinless 
hole(holon) of momentum k. 

For the application of the Holstein-Primakoff theory to functional Schrodinger pic- 
ture(FSP) representation, one has to choose the FSP operator of the charged fermion. On 
the other hand the charge-neutral fermion FSP operator is used in the slave-boson theory0. 
In the following, we derive the ground state energy and the dispersion energy of magnon for 
the two-dimensional systems of antiferromagnetically correlated electrons at and near half 
filling. 

Similarly to the case of relativistic FSP theory! we define, for chargeless scalar 
bosons(magnons) , 



> 4 -( : 



5. a) 



lb) 



with 0k, the scalar field variable in momentum space, and for charged fermions (holes), 



S 

Sul 



' & - 73 ( M k + siJ) ■ 



<5-u k J 



l.c) 



with «k, the Grassmann field variable in momentum space. The field variables and u de- 
fined above satisfy the following commutation and anticommutation relations (see Appendix 

0), 



c k ,c 



' ^k' 



J k,k', 



G?k, d\; 



(9) 



The functional Schrodinger equation is written, 



(10) 



where is the ground state for the square lattice of antiferromagnetic spin order at or 
near half-filling, 



^[0 c ,/,«] = iV^exp 



z k 



y -k"k^k 



-k"krk 



1 



k,k' 



(11) 
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where Nq is a proper normalization constant. By taking the variation of the ground state 
energy with respect to the Gaussian exponents Q c , Q d and 

^W^ 1 ^' 1 (12) 

we readily obtain Q c = Q d = 1 and Qf = — 1. For the sake of brevity, the summation symbol 
is omitted in the above expression (|TT|). 

We find from the use of Eq. (8) through (|T2| ) that the ground state energy of the total 
system is (see Appendix |A] for details) 



E = Jsz X; V 1 - Tk + 2£>k,k) - Ms + l)zN 
Jz 2 

+ 7q (-Dka.kx-q-Dkx.ka+q - -Dkx.ki-q-Dka.ka+q) (13) 

ki,k2,q 

where D k , k / = \ \(l + Q f ) + &) (l + Qf)] , 6 k , k > = (/t/k)4^- It is of note that 



k,k' 

at half-filling X) k .D k , k = 0. We readily note from Eq. ( |TB| ) above that the energy of free 
magnons at or near half-filling is given by 

E = J2^-Js(s + l)zN, (U.a) 

k 

where u>^ is the dispersion energy of magnons, 



u% = Jszy/l-^. (14.6) 
for a square lattice, the dispersion energy is simply, 



o£ = 2JV1-t£ (14.c) 

with s = |. Eq. (14. a) represents the familiar zero-point energy of the antiferromagnetic 
system^. The last term in Eq. (14. a) is the quantum contribution which is essential for 
lowering energy. 
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III. HOLSTEIN-PRIMAKOFF BOSON APPROACH OF T — J HAMILTONIAN 
BY FUNCTIONAL SCHRODINGER REPRESENTATION; 
ANTIFERROMAGNETIC SPIN POLARON 



Based on the t — J Hamiltonian, we write the antiferromagnetic spin polaron Hamiltonian 
in momentum space@, 



H AFP = E 



2s N ^ \Jsz 



CqC q + dl.dk) - V2stZ\l — 



N 



X 



1 - J- g Jj E sinh2 ) ^ E /k/k-q ^cc(k, q)c q + ^(k, q)4 



+^E(^-i- 2 /kA]-J^iv. 



H.c. 



(15) 

k \ JS * I 

Here 0;° = JszJl — y^. t is the hopping integral. 9^ is the mixing angle defined by 
tanh 29k = 7k- Qcc and qm are the coupling functions given by 



g cc (k, q) = 7 k _ q cosh 6» q - y k sinh q 
#dd(k, q) = -7k- q sinh 6> q + y k cosh 6> q 



(16.a), 
(16.b). 



The first and the last two terms in Eq. (|T5[) above are identical to the terms that appear in 
Eq.([|). It is reminded that the last two terms were not considered in Ref.(4). 

In order to properly estimate the influence of fermi-bose (holon-magnon) coupling, we 
now introduce the shifted wave functional^ associated with the scalar field, 



u] = Nq exp 



(n 



k (0 k 



klVk 



>k0 



) 2 ) + 2 M k fi k,k'«k' 



(17) 



The scalar field is shifted by a constant 4>o- With the use of Eqs.(8) and (|17|) for Eq. (|15"|), 
we obtain the ground state energy, E = E[<f>Q, 4>q] Q c , fl d , fl?}, 



e = (H AFP )/(n\n) 



E 


< 


2s 


q 




1 


+ 




1 - nj 2 ) fi q - x 



i (l - fif) - ~ (l - ^ 



+ 



'qO 



Ld 2 

; qf) 



1 - 2 L> k k - Js 2 z7V 



'■1st : t / A ( 1 _ -LI £ S inh 2 9 kl | £ (Ac,k-q + 4,k-q) 



+<7«w(k,q)-^0j o 



4s N 



+ H.c. 



k,q 



# cc (k, q)y|0 q o 



(18) 

where -D k i = | (1 + £V)(£y + f2^)~ 1 (l + £V) . We take the variation of the above ground 



k,l 



state energy with respect to the vacuum expectation value <$j, 



5E 




<50qO 




E^k 


1 - 


k 




X E ( Ac,k 


k,p 





2sN^[jsz 1 / 



0ko^k, q - V^**y 1 M - 11 ^ sinh 2 q , 



2s , \ ^ 
q \ 



-qO ~~ 



x E (^k,k- q + 5k,k-q) ^ cc (k, q)-^= 



k,q 



2s A^ , \ Js^ 
q x 



'qO 



2s^/4 fl - llE sinh2 ^c 



2s^wl fl - llE sinh ^q' 



1 



X —7= E Ac,k-q#cc(k, q) 

V2 kq 





(19) 



Thus from the above expression we readily find the vacuum expectation value, 



\fstZ\ 


/ N 


fl - 

V 4s TV 


X) q ' sinh 2 6 




< 


1 - 


1 2 y-« 

2s TV ^q' 







E jD k',k'- q 5' cc (k / ,q). 

k' 



(20) 



We note that -Dk',k'-q = for |q| ^ and fi q = fi q = 1. Realizing the vacuum to vacuum 
transition, we obtain the ground state energy from the insertion of Eq. (EH) into Eq. (|I5D, 



£ = J S2 £ M - 7 2 _ + i) zN + £ S n (k, 



(21) 



where 



E n (k,w) 



(l-ME q 'Sinh 2 fl q 



x 



fl-Uy 

|_ x 2sTV^qVJsz , 

E (^k,k- q ^c(k, q) + Dl^gUk, q)) /wj. 



(22) 
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The first two terms in Eq. (pj]) are identical to the expressions of Eq. (14. a). It is seen 
from Eqs. (pi) ) and (]22f) that the ground state energy is lowered by the self-energy of the 
holon(spinless hole) which is contributed from coupling to magnons. 



IV. SLAVE-BOSON APPROACH OF T — J HAMILTONIAN BY FUNCTIONAL 

SCHRODINGER REPRESENTATION 

Here the slave-boson functional Schrodinger representation of the t — J Hamiltonian will 
be discussed for the hole doped system of antiferromagnetically correlated electrons. In this 
section emphasis is placed on a rigorous derivation of dispersion relation and its comparison 
with the Holstein-Primakoff theory of Heisenberg Hamiltonian in the limit of half-filling. 
The t — J Hamiltonian in the slave-boson representation is given by0 

H=~tJ2 bibjflfj, - J ~ £ {flflji-.fi. ~ flfj-Jjafi-a) ~ (23.o) 
(ij)a (ij)a i<J 

Allowing a uniform hole doping rate, 5 = bity, we rewrite 

H = ~t$^2i fitrfj<r ~~ ~2 ifitrfj-irfi-ffi(T ~~ fiafj-afjvfi-a) ~ ^^^fiafia- (23.6) 

(ij)a (ij)a i(T 

Here t is the hopping strength and //, the chemical potential. bi(b\) is the annihila- 
tion(creation) operator for a spinless boson and fi a (fL)^ the annihilation(creation) operator 
for a chargeless fermion with spin a at site i. 

Following Floreanini and JackiwB, the chargeless fermion field(spinon) operator / can be 
written in terms of the Grassmann field variables u 

where a represents (i,cr). We write Gaussian functional, 

|0) = (u|*n) = N n exp(^u a Q a(3 Uf3) = N n exp^uQu (25) 

where Q is the antisymmetric 2x2 kernel matrix, and Nq, the normalization constant. For 
brevity the summation sign over the lattice sites is omitted in the above expression. We 
now introduce Eqs. (pi]) and (BS) into Eq. (24. b) to write the ground state energy, 



E 



(Q\H\Q) 



The resulting ground state energy is, (for derivation, see Appendix 0) , 



(26) 



where 



E = -t5 D aa (j, i) - ^2D aa (i, i) 

(ij)a irr 

-i { D ^(i,i)D^ a ^ a (j,j) + D aa (j,i)D- a _ a (i,j)} 



(27) 



(28) 



with a and a , either the spin up or spin down state of electron. 

We readily find from the inspection of the chemical potential terms in Eqs. (24) and 
that the number of electron rii at site i and the local magnetization m 8 at site i are given by 



Hi = D TT (i, i) + D u (i, i) = + n ib 
rrii = D n (i,i) - D u (i,i) = n l} - n^. 



(29) 



From Eq. fl29|), we obtain 



D n (i,i) = (ni+mi)/2, 
D ll(hi) = (n< — m<)/2. 

into Eq. (^7]) leads to 
E = -t5 ]T D aa (j,i) - - ]T D aa (i,i) (nj - arrij) 

(ij)a (ij)a 

-liN(l-6). 



(30) 



Insertion of Eq. 



(31) 



From the Fourier transform of expression (|3l|) above we obtain the following ground state 
energy (see Appendix |B| for verification), 



E = ±2 ]T J [(4t5f - ( Jmf] 72 + ,Pm? - //(l - S)N - 2 J(l - <5) 2 iV 



(32) 



and the dispersion relation, 
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u; k = ±2! 



^[(4M) 2 - {Jmf] 1 l + J 2 m 2 



(33) 



with 7k = [cos(k x a) + cos(k y a)} /2. 

At half-filling, i.e.,5 = 0, Eq. (£3|) leads to 



c^ k = ±2JmJl - 7 £. 



For the system of paramagnetic state, i.e., m = 0, we obtain 



co>k = ±8t57k- 



(34) 



(35) 



For the case of vanishing hopping integral or in the limit of large J, that is, tj J <C 1, we 
obtain the dispersion energy of spin waves(magnon) from Eq. (|33|), 



cj k = 2JmJl — 7 k 



Realizing from Eq. ( f29|) that m — 1 at half-filling, we find that 



(36) 



c^ = 2JJl- 7 2 



(37) 



We find from Eq. (37) above that in the limit of half-filling, that is, 5 — > 0, the slave-boson 
theory of FSP leads to the identical dispersion relation (Eq. (14.c)) of magnon obtained 
from the Holstein-Primakoff theory. 



V. SUMMARY 

In the present study we showed the functional Schrodinger representations of Holstein- 
Primakoff boson and slave-boson theories for the Heisenberg Hamiltonian and the t — J 
Hamiltonian respectively. From the use of the functional Schrodinger picture(FSP) theory 
of Holstein-Primakoff boson approach for the Heisenberg Hamiltonian we obtained both the 
zero-point energy of free magnons and the self-energy of holon(spinless holes). The FSP of 
slave-boson theory was also introduced into the t — J Hamiltonian to derive the dispersion 
relation of magnon. We find that in the limit of half-filling(<5 — > 0) the dispersion relation 
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derived from this approach is identical to the one obtained from the Holstein-Primakoff 
boson approach. Improvement over the present approach of functional Schrodinger picture 
theory is desirable to fully account for many-body effects(correlation effects) beyond the 
mean field level. 
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APPENDIX A: EVALUATIONS OF GAUSSIAN EXPONENTS, tt c , fL D AND Q F 
AND FUNCTIONAL SCHRODINGER REPRESENTATIONS OF ONE- AND 
TWO-BODY TERMS FOR BOTH BOSONS AND FERMIONS 

In real space, we define the ground state Gaussian functional of bosoni to be 

((j)\n b } = * n [0] = N n e -|//^Wn b (x- y )^(y) ^ ^ 

or in an abbreviated form, we write \I/n[</>] = e~^ nb ^. Nq is the normalization constant. 
The ground state functional in momentum space is written 

# n [0] = jVne-sM ^* (A2) 
or allowing discreteness in momentum k, 

*n[<3 = N n exp 



(A3) 



k 

where N is the total number of lattice sites. 

Introducing the boson field 0, we define the boson operators in FSP, 

c(x) = \ ( 0(x) + ' } 



y/2 V ^0( x ) 
7l( 0(X) ^, 



c t (x) = ^f0(x)--^-l (A4) 
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which satisfy the commutation relation, 



[c(x), C t(x')] = 5(x-x'). 



(A5) 



The field operators in momentum space are then 

1 1 



c(k) 
C t(k) 



1 1 

V2 (2vr) d 



dx e 



<ix e 



ikx 



X + 



-ikx 



5 



x - 



50(xW' 



(A6) 



and thus, 



5 



c t(k) = -5= u*(k) - - 



C k , c k / 



(AT) 



^k,k'- 



the use of which satisfies the commutation relation, 

Now we evaluate the vacuum state expectation of one-body term ((c kl c k2 )) in FSP. 
Using the ground state Gaussian functional 



*n[0] = N n exp 



we obtain 



50; 



<50_ k 

1 

2 



p 

-l(fi^ k + ^ k k ) tt n $. 



By using the symmetric condition, f2 6 k = f2 6 _ k , we obtain 

5 



(A8) 



Thus we find 



(0|cL lCk2 |0) - 1 (l - fi k2 ) (l + n*J / D^lM^ntf] - i (l - ^ k2 ) 4 llk2 (A9) 
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with 



/ £0 ki k2 e"^ 
5 2 



feo 5J kl 5J£ 2 
lim 



J-o 5J kl 5J£ 2 



lim 



5 



lim -fi 



'k 2 



<W 2 + ^ _1 -^^^- rnfc " lj ) 



k 2 °ki,k 2 - 



(A10) 



By inserting Eq. (|A1Q ) into Eq. ( |A9| ) , we obtain 



lic k2 |o) - \ (1 - fi k2 ) (1 + n>J iW, k2 - 1 (1 - fi k2 ) 5 kl , 



k 2 - 



For ki = k 2 , we get 



(0|4ck|0) 



-d-ni) (l + nDnl-'-Ui 



\ K 1 - « 



(All) 



(A12) 



In order to determine the value of Q b for the energy of free boson we now take the variation 
of the one-body energy term ((c k c k ) = (0|c k c k |0)) with respect to i\, 



a 



am 



(4 c k) = 7 



fi k 







(A13) 



to obtain 



n*=±i. 



Here fi k = — 1 is discarded as the ground state Gaussian functional cannot be defined. 
Finally for ki = k 2 , Eq. (\KU\) leads to 
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(0|c£c k |0) = for n b k = 1 (A14) 

as expected. In a similar manner we obtain Q k = 1 for the two-body terms. 
We note that for fermions 

f 1 ( 5 ^ 



and 



\/2 V 5ul, 



{fa, fl} — dab 



with i — a and b to represent site indices. 
We find 



(n\Q.) = J DuDv! exp[u ] \n + tt)u} = det^tt + Q). (A15) 



Realizing that 



fb\0) -> -j= (u b + -^j exp(ultt cd u d ) 

= ~^(u b + 5 bc Q cd u d ) exp(ulQ cd u d ) 
= {hd + &bd)u d exp(ulQ cd u d ) 
= (I + Vt) bd u d exp(ultt cd u d ) 

the expectation value of (P f) is 

(flu) - D ba = |<o| («; + A) (u b + _^ |n>/<n|n> 



= ^{1 + Q) ea (I + Q) bd J DuDu ] u\u d exv\u\tt + tt)u} 
1 <5 5 

= ^i 1 + ^)ea(I + ^)bdJ^-^Z [T],f]}\ v=il=0 

1 (5 5 

= -(/ + n) ea (7 + V) bd —— exp[-rj(n + H)" 1 ??] |^ =0 

= ^(/+n) ea (/ + Q) M (fi+n) e - d 1 

= ^[(z + ^^+^-^j+n)^ (M6) 
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where 



Z [r), rf\ = J DuDv) exp[ J (u^Su + rju + vfrf)] = det(S') exp (— r)S (A17) 
with S = Q + Q. Similarly we find, for the two-body term, 

(fiftfcfd) - (/ + + U) ja (i + n) df (i + n) ch (n + U)^(n + n)^( - M Pj ) 
= [(/ + 0)^(0 +n)7/(/ +n) io ][(/ + +n)^(/ +n), b ] 
-[(/ + a) df (Q + n)7/(/ + ty lb \ [(i + q) c/1 (q + + n) ja ] 

= D^D ch - D db D ca (A18) 



with i — a,b,c and d, the site index. 

With the use of Eq. (8), we obtain in momentum space, 



W = E^ 

k 

-Js 2 zN 



0J r 



2s N ^ \Jsz 



0U> + (dtd k )) +Jszj:(^--i- 2 (/i/ k > 



i---E 



2s q y Jsz 



c-1 
k 



2 (! - fi k) 



+1 (i - n£) (i + n£) ^r 1 - ^ (i - n£) + ( ^ - 1 - 2^k, k 



-Js 2 zN, 



(A19) 



{H m - rn ) E 

ki,k 2 
,0 , ,0 



a; a; \ 

I - , kl k f I ( (ct^ct^) + (4^ki4 2 rf k 2 



+2 1 + 



^kx^k 2 

(Jsz)\ 



(Jsz) 2 , 

(Ck 1 Cki)(4a d l«) 



2 ^ 

ki,k 2 



1 - 



, ,0 . ,0 ~ 



x(l-fi kl ) (l + ^J^- 1 



^(l-o(l+tt 



-- (i - fi k2 ) (i + fi k2 ) (i - n^j (fi^- 1 + s^m^s^m 
+i (i - (i + fi k2 ) 5 klik2 (i + Q ki ) n^-V* 

-iO-o ( i+n £ a )(^) a 
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+ i(l-0(l-f4,) 



+2 (l + 



i + nQ (i - n£) (i + ) n*, -1 fe" 1 + ( 



O (i + O (i - nt ) + ^ llk2 < _1 4 1)k2 ) 
(i + 5 klik2 (i + n£j n*, -1 ^* 

^ k2 ) + ki,k 2 ) 

O (i-< 



^ (i - n£) (i + n£) n£ _1 - ^ (i - n£) 



X 



(H m -h) 



ki,k 2 
2 ^ 



a; 



k 2 



C/t/ki) ((i 2 ck 2 ) + (4 2 c k2 )) 



X 



2 ^ k ^ 2 Jm2 



ki,ki 



+^(i-0(i + ^ 2 )^ 2 _1 



2 



<1 



2 iV.^ Jsz2 

ki,k 2 



1 M 51 



ki,k! 



X 



fi k2 _1 + fi C 



k 2 



+ 



nsl a - 1 + nsj a -2 



2 TV ^ ^^^i^q^^^+q^z) 

ki,k 2 

~~2~~j\[ 7q( /ki-q (^ki,k 2 +q ~~ /k 2 +q/ki) /k 2 

ki,k 2 

'~2~J^ ^2 7q ((/ki-q/k 2 )^ki,k 2 +q — (/ki-q/k 2 +q/ki/k 2 )) 
ki,k 2 

2 

Xl-^k.k - XI (-Dk2,ki-q-Dki,k 2 +q - -Dki.ki-q-Dka.ka+q) 

1, 1,- 1,_ 



" 2 jyE^q 
q 



ki,k 2 
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where 

£> k ,k' = l[(i + nf)(n + + n')] kik , = (A23) 

By noting that 

E7k = E7-k = 0, (A24) 

k k 

we obtain 

Jz 2 

Hh-h = ~YJ^ E 7q (-Dk 2 ,k 1 -q-Dk 1 ,k 2 +q ~ Ad.kj-q-Dk^ka+q) • (A25) 
ki,k 2 ,q 

By taking the variation of the ground state energy with respect to fi k (fi k ), we find that 
fi k = f2 k = 1. The resulting ground state energy is then 



E= JszJ^N^-ll + ^L) ~Js{s + l)zN 
k ^ ' 
Jz 2 

+ ~^~ Tf X/ ^q (-Dk 2 ,ki-q-Dki,k 2 +q ~ A<i,ki-q Ac 2 ,k 2 +q) (A26) 
ki,k 2 ,q 

where n k = (f^fk) — A<,k- By taking the variation of (H ) in Eq. ( |A19| ) with respect to f2 k 
we obtain f2 k = —1, which leads to .Dk,k' = at half-filling. Thus both the free hole energy 
term and the hole-hole interaction energy (the second term in Eq.( |A26| )) term vanish. 

APPENDIX B: THE DISPERSION RELATION OF MAGNONS FROM THE T-J 

HAMILTONIAN 

From Eq. (23. b) the expectation value of the t — J Hamiltonian in slave-boson represen- 
tation is 

E=(H) 

= ~ E/ (ficrfjo) ~ ~2 E/ {(fiafj—afj—cficr) ~ (fitrfj-afjcrfi—cr)) 
(ij)cr (ij)a 

-/*£<&/«*>• (Bl) 
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In the following we express both (/„/&) and (flfbftfd) above in functional Schrodinger 
picture. Using that 



we obtain 



V2 



(u b + 5 bc Q cd u d ) exp(u c Q cd u d ) 



= (5 bd + tt bd )u d exp(u c tt cd u d ) 
= (I + Q) bd u d exp(u c Q cd u d ), 



(fih) = wi«/m 

= i(J + fi) ea (J + J DuDuu e u d exp[u(Q + 
1 (5 6 

1 (5 (5 

= -(/ + U) ea (I + tybdj^jz exp[-f](n + n)- 1 7 7 ]| r?= ^ =0 

= i(/+n) ea (/ + fi) M (fi + n) e - (i 1 
= i[(/ + fi)(fi+n)- 1 (/+n)] ba 



(B2) 



where 



^oIVtV] — J DuDuexp[J (uSu + i]u + w])] = det(S) exp[—r]S 1 rj\. (B3) 



For brevity we define 



a 



ba 



(i + n) (p, + n) 1 (/ + n) 



ba 



with S = + tt, and = det^ (Q + Q). 

We express 



1 5 5 

fcfd\0) -> -(w c + ^)K* + exp(w e fi e/ M e ) 

1 <5 

= 2^ c+ + tydfUfexp(u g Q gh u h ) 



(B4) 



19 



(I + Q) d f{u c Uf + -jj-r-Uf) axp(u g Q gh u h ) 
(I + Q) d f(u c u f - UfQchii^e^ 

(I + Q) df (I + Q) ch u h u f e Wu , (B5) 



and 



(o\fiflf c f d \o) = \(i + ty lb (i + U) ja (i + n) df (i + n) ch 



x J DuDue mu u j u l u h u f e mu (B6) 



where 



5r}j 8r]i 5r) h 8r]f 

fymnVn] \rj=fi-- 

8 & & 

Art 1 ( ^ I U)8 flll ( Q + U)^ n — — —r] n exp[-rf (n + H)" 1 ^] | ^=^=0 



1 (5 (5 & & 

-det^^ + n)— -— — expH- m (fi+n)-^]|^ =0 



1 (5 (5 

= -deta(fi + 0)5 /n A o (0 + 0)~* (O + ^l)~pj^-^-r] n r] p exp[-r] q ( y n + f2)~V]|r,=rj=o 

= -det^n + U)(n + 0)^(0 + U)-^(5 nlVp - 5 plVn )e-^ a+ ^\^ 

= _det^(fi + + fi^O + Tl) h p(5 n i5 pj - S pl 8 nj ) (B7) 
Realizing that = det^(fi + we find from Eq. (|B|) and Eq. flB7|) that 



(flflfcfd) = (I + tylb(I + H)ia(/ + 0)^(1 + Q) c/l (Q + H^ft + ^^(V^ - Mpi) 

= [(/ + n) df (n + 0)7/(1 + U) ja ] [(/ + n) ch (n + n)^(/ + o) Z6 ] 
-[(/ + n) df (n + n)7/(i + U) lb ] [(i + n) ch (n + n)^(i + H) ia ] 

= D da D cb — D db D ca . (B8) 



The substitution of Eqs. (|B^ ) and (|B8| ) into Eq. leads to 



E = -t5J2D™tiJ) ~ ^E - D^ <7 (j,i)D^ a (i,j)\ 

(ij)a (ij)a 

- [D_ w (i,i)D <T _ ff C7,j) - D aa (j,i)D^ a (i, j))} - ^D a(T (i,i). (B9) 
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Allowing the global SU(2) symmetry and thus 

D a _ a (i,j) = (flj l(T ) = 0, (BIO) 

Eq. (p9| ) is rewritten, 

E = -t5^2D a(T (j,i) {D m {i,i)D_ a _ a (jJ)+D m {j,i)D- ir - <r (i,j)} 

(ij)a (ij)a 

-fT£D**(i,i). (Bll) 

From the inspection of the last term, we note that the number of electron rij at site i 
and the local (site) magnetization rrii at site i are given by 

n { = £>TT(M) +-°U(M) = n i\ + n ili ( B12 ) 

m t = D TT (i,i) -D u (i,i) = n A - n n . (B13) 



From Eqs. ( pT2|) and (|BT^ ), we obtain 

£> U (M) = -(rii-mi), 
J2 D aa(i,i) = D^(i,i)+D u (i,i) = th 
J2D aa (i 1 i)D_ a _ a (j,j) = Dtf(i,i)D u (j,j) + D u (i,i)D^(j,j) 

(7 

s (tia — rrin) _ , . (n,- + 



-^2D aa (i,i)(nj - arrij) (B14) 



with cr = 1(— 1) for up (down) spin. Using Eq. ( B14|) , we rewrite Eq. ( Bll|) 



(ij)a (ij) (ij)a 

(B15) 

i 

Allowing uniform hole doping and thus n = rii = 1 — 5, we obtain the total number of 
electrons 
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$> = XX 1 - *) = tfC 1 - *) 



and 



E 



53 



(ij) (*i> 

With the use of Eqs. ( jBT§ and QET7D , Eq. JBI5| ) leads to 



(ij) 

-JN(l-5) 2 - fiN{l-S). 



J 



Dam 



We define 



Then we have 



= D^(i,j) + D u (i,j), 
m(ij) = D^(i,j) - D u (i,j). 



D u(hj) = \[n(i, j) - m(i, j)}. 



Using Eq. (|B20|) , we write the third term in Eq. (P18| ), 



^ D a<r {j,i)D. a - (i,j) =E[ £) nO'>i)^u(i ) i)+^uO'>i)^tt(*^*)] 

foV (ij) 

{n(i,j) - m(i,j)} 



E 

(ij) 



{n(i,j) + m(i,j)} 



E D ccr{j,i)[n{i,j) -am(i,j)], 



(ij)o- 



The substitution of Eq. (g2Tj) into Eq. flBlq) leads to 

£ = E D ^(j,i) + j E D aa (i,i)am j - 7 E D ^(j,i) [n(ij) - am(i, 
-JN(1 -5) 2 -fiN(l -5). 
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We note that the two-dimensional staggered magnetization is represented by mi 
with Q = (ft, ft) and i = (i x ,i y ). Thus Eq. (|B22|) leads to 



Jul 

E = -tS E D aa (j,i) + — E oD aa (i,i 

(ij)a (ij)a 



)e iQ - i -^D aa (j,i)[n(iJ) 



am 



-fiN(l - 6) - JN(1 - 6)' 



where 



Introducing 



and 



J2 D <ra(j,i) = D n (j,i) + D u (j,i) = n(j,i), 

cr 

E D <ra(j,i)n(i,j) = ^n(j,i)n(i,j). 

feV (ij) 



E D aa (j,t)n(t,j) ~ E^ = 4Ar (! " S ) 2 

(ij) 



m(i,j) ~ rrii = me 1 ®' 1 , 



we obtain from Eq. (|B23 ) 



E = -t5j2 Daa(j,i) + ^ E crD aa (i,i)e^ + ^ E oD„&i)e«** 

(ij)a (ij)a (ij)cr 



-fiN(l-5)-2JN(l-5) 2 
-HN(1 - 5) - 2JN(1 - 5f. 



+ ^aD a Ji,i)e^ 



In momentum space, we note that for the case of square lattice, 



D„& = (23a E <U k ' k')e ik j e - ik, i , 



k,k' 



and thus 



J2D aa (j,i) 

(ij) 



(2vr) 2 



(2vr) 2 



<ii> k,k' 

E^^k^OE^e-*'' 1 . 
k,k' (ij) 
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Using 

^ e ik-j e -ik'.i = ^ e ik.(i+a) e -ik'.i 
(ij) i,a 

= E e i(k_k ' ) ' i E e ika 

i a 

= (2vr) 2 5(k - k') x 2(cos k x a + cos k y a), (B30) 
with a, the four nearest neighbour sites around site i, we obtain from Eq. 



E D aa (j, i) = E daaiK k')8(k - k') x 2(cos k x a + cos k y a) 
(ij) k,k' 

— d aa (ia, k)2(cos k x a + cos k y a) (B31) 



Noting that 



E £>„(,', z)^ = E E <U k , ky^'V^ 



(y) v-v (y)k.k' 



E <U k , k ') E e^^e^'V^ 

\ Z7Y ) k,k' i,a 

1 E^^kOE^^'E^ 

k,k' i a 



(2vr) 2 



E d aa (k, k')5(k — k' + Q)2(cos k x a + cos fc. 



a 



k,k' 



E ^0-0- (k, k + Q)2(cos k x a + cos k y a), (B32) 



and 



(ij) v / i,a k,k' 



l — 2 E d aa (k, k') e e^-^- 1 E e iQ ' E 

^ k,k' i a 



(2 

= -4E^(k,k+Q),) (B33) 

k 

we obtain, with the use of Eq. ( |B31| ) above, 

Jul 

E = —tS E 2(cos k x a + cos k y a)d aa (k, k) H E c(cos /c x a + cos k y a)d aiJ (k, k + Q) 

k,cr k,cr 

-Jm E ^( k , k + Q) - fiN(l -5)- 2JN{1 - 5) 2 . (B34) 

k,CT 
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Introducing 7k = (cos k x a + cos k y a) /2, the above equation is rewritten 
£ = E [-4t57kGU(k, k) + Jma ( 7k - 1) d aa (k, k + Q)] 



k,cr 



-fiN{l - 5) - 2JN{1 - S) 2 . 



(B35) 



We take k -sum in the reduced Brillouin zone i.e., half the first Brillouin zone, to write 
E 7iA CT (k, k) = E [7krfaa(k, k) + 7 k+Q d ff(J (k + Q, k + Q)] 

(B36) 



k.c 



k.c 



7k E [<M k > k) - cU(k + Q, k + Q)] 

k,cr 



and 



E (7k - 1) Cx(k, k + Q) = E [(7k - 1) cU(k, k + Q) + ( 7k+Q - 1) d aa (k + Q, k + 2Q)] 

k,cr k,cr 

= E t(7k - 1) cUk, k + Q) - ( 7k + 1) <C(k + Q, k)] . (B37) 

k,(T 

where we used the nesting condition, 7k +Q = — 7k . The symbol ' indicates momentum 
summation in the reduced Brillouine zone. By applying ( B36|) and ( |B37|) to (|B35| ), we 
obtain the ground state energy of the two-dimensional hole-doped systems, 



k,<r 



E = E [-4t5 7 kGU(k, k) + Jma ( 7k - 1) d aa (k, k + Q) 

-Jma ( 7k + 1) d aa (k + Q, k) + 4t5 7k cU(k + Q, k + Q)] 
-fiN(l - 6) - 2JN{1 - 5f. 



(B38) 



Realizing the equivalence between d aa (k, k') and / k // k ' from the inspection of Eq. (23. b) 
and Eq. fl2"7|) , we consider the transformation matrix in 



ik /k 



k+Q 



-4t5 7k Jmcr ( 7k - T 
-Jma ( 7k + 1) 4t5 7k 

From the determinant of the above matrix, we obtain 



ft 



k+Q 



e i -(4^7 k )' + J 2 m 2 ( 7 ^-l)=0. 



The dispersion energy is then 



(B39) 
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= ±v/(4t(57 k ) 2 - J 2 m2 (72-I). 



(B40) 



Considering summation over the up and down spins, we write the dispersion energy, 



= E v[( 4 ^) 2 - ( Jm ) 2 } 7k + J 2 ™ 2 



2 J[(U5) 2 - (Jm) 2 ] 72 + ,Pm 2 . 



(B41) 



The ground state energy of the hole-doped two-dimensional antiferromagnet is then 



E k = ±2 ^ J [(U5) 2 - (Jm) 2 ] 72 + .Pm 2 - -5)- 2 J(l - 5) 2 



(B42) 



For the paramagnetic states, m = 0, we obtain 



u; k = ±4t<57 k . 



(B43) 



For the undoped antiferromagnet, that is, 5 = 0, the dispersion energy of the antiferro mag- 
netic magnon is 



= 2Jm^fl 



(B44) 
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